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ABSTRACT 

The sub-mm bump observed in the spectrum of Sgr A* appears to indicate 
the existence of a compact emitting component within several Schwarzschild 
radii, r s , of the nucleus at the Galactic Center. This is interesting in view 
of the predicted circularized flow within ~ 5 — 10 rs, based on detailed 
multi-dimensional hydrodynamic simulations of Bondi-Hoyle accretion onto 
this unusual object. In this paper, we examine the physics of magnetic field 
generation by a Keplerian dynamo subject to the conditions pertaining to Sgr 
A*, and show that the sub-mm bump can be produced by thermal synchrotron 
emission in this inner region. This spectral feature may therefore be taken 
as indirect evidence for the existence of this circularization. In addition, the 
self-Comptonization of the sub-mm bump appears to produce an X-ray flux 
exceeding that due to bremsstrahlung from this region, which may account for 
the X-ray counterpart to Sgr A* discovered recently by Chandra. However, the 
required accretion rate in the Keplerian flow is orders of magnitude smaller 
than that predicted by the Bondi-Hoyle simulations. We speculate that rapid 
evaporation, in the form of a wind, may ensue from the heating associated with 
turbulent mixing of gas elements with large eccentricity as they settle down into 
a more or less circular (i.e., low eccentricity) trajectory. The spectrum of Sgr A* 
longward of ~ 1 — 2 mm may be generated outside of the Keplerian flow, where 
the gas is making a transition from a quasi-spherical infall into a circularized 
pattern. 

Subject headings: accretion — black hole physics — hydrodynamics — Galaxy: 
center — magnetic fields: dynamo — magnetohydrodynamics — plasmas 
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1. Introduction 

The non-luminous matter concentrated within the inner 0.015 pc at the Galactic 
Center has a measured mass of 2.6 ± 0.2 x 10 6 M (Genzel et al. 1996; Eckart & Genzel 
1996; Eckart & Genzel 1997; Ghez et al. 1998). Most of it appears to be associated with Sgr 
A*, a bright, compact radio source (Balick & Brown 1974) that anchors the stars and gas 
locked in its vicinity, and provides possibly the most compelling evidence for the existence 
of supermassive black holes. 

The spectrum of this unusual object can be described as a power-law with an index 
a that varies within the range 0.19 — 0.34 (S u oc is a ) from cm to mm wavelengths. 
However, one of the most interesting features currently under focus is the suggestion of a 
sub- millimeter (sub-mm) spectral bump (Zylka et al. 1992; Zylka et al. 1995), since the 
highest frequencies appear to correspond to the smallest spatial scales (e.g., Melia, Jokipii 
& Narayanan 1992; Melia 1994; Coker and Melia 2000). In the case of Sgr A* one expects 
the sub-mm emission to come directly from the vicinity of the black hole. The existence of 
this bump (or "excess") has been uncertain due to the variability of Sgr A*, but is now well 
established following a set of simultaneous observations (from A20cm to A 1mm) using the 
VLA, BIMA, Nobeyama 45 m, & IRAM 30 m telescopes (Falcke, et al. 1998). 

The behavior of Sgr A* is dictated by the manner with which plasma accretes onto 
it from the nearby environment. Three-dimensional hydrodynamic simulations of the gas 
dynamics at the location of Sgr A* (Coker & Melia 1997) indicate that the accreted specific 
angular momentum A (in units of crs, where r$ = 2GM/c 2 is the Schwarzschild radius 
in terms of the black hole mass M) can vary by 50% over ^ 200 years with an average 
equilibrium value for A of 40 ± 10. Thus, even with a possibly large amount of angular 
momentum present in the wind surrounding the nucleus, relatively little specific angular 
momentum is accreted. This is understandable since clumps of gas with a high specific 
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angular momentum do not penetrate to within the capture radius, Ra = 2GM/v w 2 , defined 
in terms of the wind velocity v w at infinity. The variability in the sign of the components of 
A suggests that if an accretion disk forms at all, it dissolves, and reforms (perhaps) with an 
opposite sense of spin on a time scale of ~ 100 years. 

The captured gas is highly ionized and magnetized, so it radiates via bremsstrahlung, 
cyclo-synchrotron and inverse Compton processes. However, this emissivity appears to be 
inefficient in the case of Sgr A*, so most of the dissipated energy within the large scale 
quasi-spherical inflow is carried inwards (Shapiro 1973; Ipser & Price 1977; Melia 1992). 
A viable explanation for Sgr A*'s low radiating efficiency is that the advected magnetic 
field is well below its equipartition value. This may not be surprising in view of the fact 
that the actual value of B depends strongly on the mechanism of field line annihilation, 
which is poorly understood. Two processes that have been proposed are (i) the Petschek 
(1964) mechanism, in which dissipation of the sheared magnetic field occurs in the form 
of shock waves surrounding special neutral points in the current sheets and thus, nearly 
all the dissipated magnetic energy is converted into the magnetic energy carried by the 
emergent shocks; and (ii) van Hoven's (1979) tearing mode instability, which relies on 
resistive diffusion of the magnetic field and is very sensitive to the physical state of the gas. 
In either case, the magnetic field dissipation rate is a strong function of the gas temperature 
and density, so that assuming a fixed ratio of the magnetic field to its equipartition value 
may not be appropriate. 

Kowalenko & Melia (1999) have used the van Hoven prescription to calculate the 
magnetic field annihilation rate in a cube of ionized gas being compressed at a rate 
commensurate with that expected for free-fall velocity onto the nucleus at the Galactic 
Center. Whereas the rate of increase dB/dt\f in B due to flux conservation depends only 
on the rate f of the gas, the dissipation rate dB/dt\d is a function of the state variables and 
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it is therefore not necessarily correlated with r. Although these attempts at developing a 
physical model for magnetic field dissipation in converging flows is still rather simplistic, it 
is apparent from the test simulations that the equipartition assumption is not always a good 
approximation to the actual state of a magnetohydrodynamic flow, and very importantly, 
that the violation of equipartition can vary in degree from large to small radii, in either 
direction. 

The first serious attempt at modeling the spectrum of Sgr A* using a sub-equipartition 
profile for the magnetic field was made by Coker & Melia (2000), who adopted a spherical 
infall as a simplified version of the actual accretion picture. Of course, the real accretion 
flow will deviate from radial at small distances from the black hole, where the gas begins 
to circularize with its advected specific angular momentum. This deviation of the accreting 
gas away from a purely radial infall may provide a clue for the appearance of the bump in 
the spectrum at sub-mm wavelengths, which seems to hint at the existence of a distinct 
geometry for the emitter at this energy. We here suggest that the sub-mm "excess" in the 
spectrum of Sgr A* may be the first indirect evidence for the anticipated circularization of 
the gas falling into the black hole at 5 — 25 r s . Although the physical conditions within 
the quasi-spherical infall evidently suppress the magnetic field well below its equipartition 
value, this need not be the case once the gas circularizes and forms a Keplerian structure. 
Indeed, it is expected that a magnetic dynamo within the differentially rotating region may 
overwhelm the field annihilation rate and actually lead to a saturated field intensity. 

It is our intention here to fully explore the magnetic properties of this inner circulating 
region, and to assess the likelihood of producing the sub-mm spectral bump with a magnetic 
dynamo. Given the wide latitude of possible configurations in the outer quasi-spherical 
infall for producing the longer wavelength radio emission, we will not here attempt to 
construct a complete picture for the whole accretion region. Rather, this task is better 
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coupled to a hydro dynamic simulation that self-consistently merges the large scale flow to 
the circularized structure at smaller radii. These calculations are now in progress, and their 
results will be reported in the future. 

We will first summarize the key physical principles underlying the dynamo process, 
and then examine what configuration the magnetic field should have in the Keplerian 
region surrounding Sgr A*. We will then demonstrate that an excellent fit to the sub-mm 
data is possible with this picture, and discuss the implications of this model to our overall 
understanding of the environment surrounding the massive black hole at the Galactic 
Center. 

2. The Magnetohydrodynamic Dynamo in a Rotational System 

Since the discovery of the local shear instability in weakly magnetized disks (Balbus 
& Hawley 1991, hereafter BH1), several numerical simulations (Hawley & Balbus 1991, 
hereafter HB1; Hawley & Balbus 1992, hereafter HB2; Hawley, Gammie & Balbus 1995, 
hereafter HGB1; Stone, Hawley, Gammie & Balbus 1996, hereafter SHGB) have confirmed 
the fact that this instability plays a crucial role in rotational accretion systems. A linear 
analysis (Balbus & Hawley 1992, hereafter BH2) has demonstrated that the instability is 
extremely powerful. Its maximal growth rate is of the order of the angular velocity, provided 
that the latter decreases outward and that initially there exists a weak magnetic field. 
In the case of an axisymmetric perturbation with a weak B z component of the magnetic 
field (BH1), the maximal growth rate in a Keplerian disk reaches 0.75 Q at k z ~ Q/vaz, 
where Q is the angular velocity, k z is the ^-component of the perturbation wavenumber 
and vaz = \f^ll 1S the Alfven speed in the z-direction. The numerical simulations 
show that the instability saturates at a turbulent state, producing a significant angular 
momentum flux, which is dominated by the Maxwell stress rather than the Reynolds stress. 
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This process has been invoked to account for the origin of the anomalous viscosity in these 
systems (SHGB). 

Another important consequence of this instability, one that we explore at length 
in this paper, is the amplification of the magnetic field. It has been shown that this 
instability constitutes a magnetohydrodynamic dynamo (Hawley, Gammie & Balbus 1996, 
HGB2; Brandenburg, Nordlund, Stein & Torkelsson 1995, hereafter BNST). An external 
magnetic field, which is sometimes invoked to magnetize the disk, is not necessary for the 
instability to work. An internal turbulent magnetic field can also drive the instability. 
Once one introduces the complete magnetohydrodynamic equations to describe a shearing 
hydromagnetic structure, the system saturates at a turbulent state that includes a significant 
magnetic field energy density. Numerical difficulties with the simulations have thus far 
prevented the acquisition of quantitative results that may be used directly to describe a 
real astrophysical situation. However, based on what is now understood about the dynamo 
process, some very important qualitative results have been obtained, and these can provide 
a guide to the manner in which the magnetic field grows within the converging flow around 
a compact object. 

2.1. The Basic Equations 

For the sake of completeness, let us first analyze how the instability develops in a 
weakly magnetized Keplerian flow. The basic dynamical equations are (BH1): 



d\np 



+ V- v 



(1) 



dt 






dB 
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where d/dt is the Lagrangian derivative and $ is the external gravitational potential. The 
other symbols have their usual meaning. We will adopt standard cylindrical coordinates 
(r, 0, z), where r is the perpendicular distance from the z-axis. 

Because the maximal growth rate is reached in the axisymmetric case with a weak 
vertical field, i.e., with B = (0, 0, B z ), we will here consider perturbations for this specific 
situation. As such, the Eulerian perturbations, which we denote by 5v, 5B etc., are 
modulated by the function e »(*vr+fc*z-<jt) ^ w ] iere an( j are; respectively, the radial and 
vertical components of the wavevector. The numerical simulations show that buoyancy is 
not a significant factor influencing the instability, nor is the compressibility of the fluid. By 
neglecting these terms and assuming incompressibility, so that 5p = in all the equations 
other than the equation of motion and the equation of state (Balbus & Hawley 1991), 
and that 5P = in the equation of state, one obtains the following linearized dynamical 
equations: 

k r 5v r + k z 5v z = 0, (4) 

dt Anp \ p Anp ) 

^ - -ik 5 ^- (6) 
dt ~ %kz p ' {b) 

-dT + 2n 6Vr = % ^T P bB ^ (7) 

d8B r 



dt 
dSB z 

dt 



ik z B z 5v r , (8) 
ik z B z 5v z , (9) 
5B r + ik z B z 5v (f> , (10) 



dt dr 

where Q is the angular velocity in the circularized flow, and k 2 = (2Vt/r) d(r 2 fl)/dr is the 
square of the epicyclic frequency. 

Replacing the Lagrangian derivatives with respect to time t by — iu in the linearized 
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equations and eliminating the Eulerian perturbations, one obtains the dispersion relation 

~ k 2 z v\ z f - % «V - kl v 2 Az ) - 4fi 2 M^i = o, (11) 



where k 2 = k 2 + k 2 . In the case of Keplerian flows, k — tt. This equation can easily be 
solved for u, which yields 



, .2 _ J 2 , K z0 nJ z0 I z0 (12) 

UJ o-k z0 + 2k , 2y k , + 16k y W 

where u = u/Q and k and k z0 are, respectively, k and k z expressed in units of Q/va z - It is 
noted that u 2 reaches its minimum value of — 9Q 2 /16 when k 2 = k 2 = (15 / 16) (Q/va z ) 2 ■ For 
k r = 0, the modes become stable when k 2 > 3 (Q/vaz) 2 , and in the long wavelength limit 

cu 2 ~-3(v Az k z ) 2 . 

In order to appreciate the physical implication of this instability, we examine the 
fastest growing mode, which occurs when k 2 = k 2 = (15/16)(Q/v a z ) 2 an d uj 2 = — 9Q 2 /16. 
Solving the linearized equations (p|)-(|TUD, we get 

5v r = 5v<f,, (13) 
5B r = -5B+, (14) 
SB r = i^-^f-for, (15) 

sb, = -4|f (-) 



\5B r \ 2 5p\5v r \ 2 



8tt 3 2 

Keeping the time derivative terms in the linearized equations (|4])-(|i"0D and using this 
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solution, we interpret the physics of this unstable mode as follows: The perturbation 5v r , 
which is generated by 5v<j, through the Coriolis force term 2Q5v ( j ) in Equation (||), induces 
the perturbation 5B r through Equation (|||). Due to the shearing of the disk, the term 
r dQ/dr 5B r in Equation (fl0|) then shows that SB r leads to the production of a perturbation 
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SB^, which in turn enhances Sv^ through the right hand side of Equation (0). Thus a 
positive feedback loop is established. Some of the other terms in the linearized equations 
act to stabilize the perturbation, but they are overwhelmed by the positive feedback in the 
unstable modes. However, for the modes with a large wavenumber, the term ik z B z bv^ in 
Equation ( |10D and the term ik z B z 5B r / Anp in Equation @ will overwhelm the positive 
feedback and make the mode stable. 



2.2. The Dynamo in Sgr A* 



The solution represented by Equation fll7Q deserves special attention; it says that the 
turbulent kinetic energy density is approximately equal to the turbulent magnetic field 
energy density. That is, it points to an equipartition of kinetic and magnetic field energy 
densities in the final saturated state of the system. Several numerical simulations (BNST, 
HGB1, HGB2, SHGB) have largely confirmed this result. They show that the ratio of these 
energy densities in the final turbulent state of the system is only weakly dependent on the 
initial and subsequent physical conditions assumed in the calculations. We shall therefore 
here adopt the conclusion drawn from these studies, that in rotational flows 

<5B 2 > 1 



C - < P 5v 2 > , (18) 

where the constant Cq has a value between 1 and 10, depending on the vertical profile of 
the Keplerian structure. 

Although the above analysis is based on a very specific model, many numerical 
simulations (BNST, HGB2) have demonstrated that this instability exists more generally, 
and that even without an external magnetic field, it can produce a significant magnetic 
energy density, which therefore constitutes a hydromagnetic dynamo. We are here primarily 
interested in the final saturated state, whose existence has been shown to be inevitable by 
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the numerical simulations. In the following, we will build a simple picture for the dynamo 
process in Sgr A*, and show how the turbulence is sustained. From these simulations, 
we infer that the vertical component of the internal turbulent magnetic field can produce 
unstable wave modes similar to those produced by an external magnetic field and that the 
dynamo is driven by them. Since we adopt a zero value for the external magnetic field 
(a reasonable approximation consistent with the view that the turbulent magnetic field 
exterior to the Keplerian flow is greatly sub-equipartition, Kowalenko & Melia 1999), the 
turbulent magnetic field is in fact then the total field B. Applying the above analysis to 
the situation in Sgr A*, we replace the external vertical field with the mean-square root 
( a/ < 5B 2 > = a/ < B 2 >) of the turbulent magnetic field's vertical component (hereafter 
the symbol <> denotes the mean value of the particular physical quantity). However, 
because the final state is turbulent, it is expected that the growth rate will be smaller than 
that in the case where there exists an underlying ordered external field. 

Before proceeding with the model, we need to examine whether or not this instability 
can be damped by Ohmic diffusion for the conditions expected in Sgr A*. Jin (1996) first 
discussed the possible damping of the shear instability using linear analysis. More recently, 
this issue was addressed numerically by Fleming, Stone and Hawley (2000). The overall 
conclusion is that the instability is effectively damped when the diffusion length becomes 
comparable to the wavelength of the most unstable mode within the latter's growth time 
scale, i.e., roughly one revolution period. As we have seen, the most unstable wave mode 
in the Keplerian flow is k 2 = k 2 z = lhVt/lQv Az , where va z = \/<B\> /An < p > is the 
vertical Alfven speed of the turbulent magnetic field. As the magnetic field is relatively 
weak, vaz is small, and all the modes with k z < 3Q/va z are unstable. These unstable modes 
lead to an increase in < B 2 Z >, which in turn leads to an increase in the value of va z - So the 
wavelength \ z = 2ir/k z of the most unstable mode increases. There does not appear to be 
any physical reason why the magnetic field should stop increasing before \ z has reached a 
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value equal to the height H of the Keplerian flow. So the length scale of the most unstable 
wave mode is expected to be H. 

Let us now show that the diffusion length scale within one period in Sgr A* is much 
smaller than this wavelength. Using the parameter values characterizing the accretion onto 
this object, we find that the angular velocity of a Keplerian flow near ten Schwarzschild 
radii (r s = 2GM/c 2 , where for Sgr A*, M « 2.6 x 10 6 M ) is « ^GM/{lQr s f ^ 0.001 
s -1 . The temperature of the plasma in this region is expected to be around a few 10 10 K 
(Melia 1994; Coker & Melia 2000), so this gas is fully ionized, and its resistivity is (Lyman 
1961) 

r] = 7.26 x 10" 9 ^ (c.g.s.) , (19) 
T2 

where the Coulomb logarithm In A is a slowly varying function of the electron density and 
temperature; 30 is a reasonable value for the conditions at the Galactic Center. Thus, for 
these parameter values, rj ~ 0.22 x 10 -21 (c.g.s.). The diffusion length scale corresponding 
to this resistivity is 

'-fig- < 20 » 

where r is the diffusion time, which we set equal to the growth time scale of the most 
unstable wave mode, i.e., about 1000 s here. Thus, we infer a diffusion length of about 
4 cm, which is clearly much smaller than any reasonable value of H, and therefore the 
characteristic wavelengths of the unstable modes. We conclude that the dynamo in Sgr A* 
cannot be damped by Ohmic diffusion. 

One of the important properties of the final turbulent magnetic field is that it is 
dominated by its azimuthal component. All the simulations show that the azimuthal 
component counts for about 80 percent of the total magnetic field energy. This result can 
be understood by examining the equations describing the evolution of the magnetic energy 
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density (remembering still that the turbulent field generated by the dynamo constitutes the 
total field in this system): 

YmT = r ^^ + ^^ x ( 5vxB )U + ^^l v2B U' ( 21 ) 

1 BB 2 - nr 2 

— = BJV x (5v x B)] r + J— 5JV 2 BL , (22) 

2 dt L v ;J 4tt 1 1 v ; 

1 r)ft 2 nr 2 

= B z [Vx(5vxB)] z + ^-B z \V 2 B\ z , (23) 

2 Ot 47T 

where is again the resistivity of the plasma, and c is the speed of light. A linear analysis 
shows that the amplitudes of the azimuthal and radial components of the magnetic field 
are equal when the perturbation is small, but that the final turbulent state is affected by 
the nonlinear character of the magnetohydrodynamic equations. As the amplitudes of the 
perturbation increase, nonlinear effects become more important. Due to the shearing of 
the Keplerian flow, the average value of rB^B r d£lj dr is positive. The energy Equations 
(HI) — (HI) therefore show that rB^B r d£lj dr contributes to a growing anisotropy of the 
turbulent magnetic field, in the sense that more and more magnetic field energy is generated 
in the azimuthal direction. For a Keplerian flow with rdQ/dr = —3Q/2, the growth rate 
due to the shearing of this structure is larger than that associated with any other dynamo 
process. So the azimuthal component of the magnetic field dominates the final magnetic 
field energy density. 

It is interesting to note that this simple scenario may also be used to interpret 
the anisotropic spectrum of the turbulent magnetic field energy density observed in the 
numerical simulations (HGB1). Because the radial magnetic field is stretched in the 
azimuthal direction, the magnetic energy carried by the modes with the largest wavelengths 
should be associated mostly with their azimuthal components. However the spectrum for 
the vertical and radial components should be similar to each other. 

So the situation with regard to Sgr A* is the following: as the gas flows inwards 
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and spirals into an approximately Keplerian structure at a distance from the black hole 
corresponding to the circularization radius, a linear instability first stretches the magnetic 
field lines carried by the gas and produces a radial component of B. The magnetic field 
generated during this step is approximately in equipartition with the turbulent kinetic 
energy and counts for a small fraction of its final intensity, but it nonetheless provides 
the seed for the next step. Second, the shearing in the Keplerian flow stretches the radial 
magnetic field in the azimuthal direction and it increases the magnetic field energy (as seen, 
e.g., in HB2). The energy comes from the rotational energy of the gas, and during this 
process a significant amount of angular momentum can be transported outward as a result 
of the Maxwell stress BrB^/Air. Finally, some of the magnetic field energy is converted 
into kinetic energy, which is eventually dissipated into thermal energy as a result of the 
viscosity, through the Lorentz force term B ■ [V x (<5v x B)]. In addition, the magnetic field 
energy may be dissipated through Ohmic resistivity. For example, with their assumed large 
numerical resistivity, BNST found that about half of the magnetic field energy is dissipated 
in this way. However, in the case of the Galactic Center, the actual resistivity is rather 
small and rjc 2 /Awrs — 10~ 13 cm s" 1 is insignificant compared to any velocity scale within 
the plasma. The dissipation of the magnetic field in this fashion is therefore negligible. 

So from the magnetic energy Equations fl2"ip — (P3"[), we infer the following proportionality 
based on dimensional analysis: 



which is the main result we have been seeking. Note that for a Keplerian flow, 




(24) 



where H is the height of the flow and we have used \J < p 5v 2 > / y/< p > to characterize 
the turbulent velocity. Combining Equations fll8|) and (|24]) , we therefore obtain 




dr 



Q, = (GM/r 3 ) 1/2 , for which rdtl/dr 



(3/2) a 



- 15 - 



Some support for the validity of this relationship has been provided by the simulations 
reported in BNST, which include an analysis of the magnetohydrodynamic dynamo at two 
different radii, one of which is five times smaller than the other. Since this Keplerian flow 
is characterized by the following equations, 

(Hn) 2 oc P/p, (26) 
p oc r 15/8 , (27) 
H oc r 9 / 8 , (28) 

the expression in ( P5|) suggests that < B 2 >oc P, the gas pressure. An inspection of the 
numerical results in BNST shows that the ratio of magnetic field energy density to thermal 
energy density at the two radii is 0.013 and 0.014. Given that the ratios of the densities, 
the heights and the angular velocities at the two radii are, respectively, 25, 5, and 10, 
the constancy of the ratio of magnetic field energy density to thermal energy density is a 
strong indication that Equation ( ^5|) provides an adequate representation of the saturated 
magnetic field intensity under these conditions. The validity of this expression is also 
supported by the global simulations of Hawley (2000). 

We conclude that the magnetic field intensity within the gas converging onto Sgr A* 
approaches the functional form given in Equation (P5|) once the flow settles into a Keplerian 
structure at small radii. The fluctuations associated with the accreted angular momentum 
(Coker & Melia 1997) correspond to a circularization radius ~ 5 — 50 r$, so we anticipate 
that B will approach the distribution in Equation fl25|) within this region. 



3. Calculation of the Spectrum 
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3.1. The Anomalous Viscosity 

In a Keplerian flow with column density £ and angular velocity f2 = (GM/r 3 ) 1 ' 2 , the 
radial velocity v r at (cylindrical) radius r is given as (e.g., Stoeger 1980) 



where v is the kinematic viscosity, 

2 W r4> 

and W r( f, is the vertically integrated sum of the Maxwell and Reynolds stresses (Balbus et 
al. 1994). For the problem at hand, the Maxwell stress dominates, and 



(the average inside the integral being taken over time). Even though this approximation is 
valid simply on the basis that the Reynolds stress is relatively small, its validity is enhanced 
by the fact that since the turbulent velocity (which accounts for this kinetic stress) and B 
are generated by the same process, both should be scalable by B. Numerical simulations 
(e.g., by Brandenburg, et al. 1995) show that (3 V changes very slowly with r. In the 
particular cases considered by these authors, (3 U ranged in value from m 0.1 to 0.2, which 
represents an increase by a factor of only 2 as r decreased by a factor of 5. We will here 
find that a "mean" value of ~ 0.2 — 0.3 is required for this quantity. 

To use the result of § 2.2, we first need to know the vertical structure of the Keplerian 
flow. For steady conditions, one can obtain the vertical profile by assuming that the gas is 
(on average) in local hydrostatic equilibrium. Balancing gravity and the pressure gradient 
in the vertical direction, we obtain the scale height 



3 d 
r 1 / 2 S dr 




1/2 



(29) 




(31) 




(32) 
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where T is the gas temperature at radius r, R g is the gas constant, and /i is the molecular 
weight. As an approximation, we will assume that the Keplerian flow is axisymmetric and 
is independent of the vertical coordinate. Written another way, we have 

(Htt) 2 = 2P/p . (33) 

Following the result of § 2.2, we will write 

/d2 r p y tr p 

^ w pdz =Pp^jr> ( 34 ) 

where (3 P is roughly constant with a value of ~ 0.03. In our simulations, we will find that a 
value ~ 0.03 — 0.04 is required to match the observations. Thus, with 

M = — 27rr£tv , (35) 

we can integrate Equation ( p9[ ) to obtain 

v r = rT ^\ /9 , (36) 

T r /v r0 + ^(GM)V2( r j/2 _ r i/2y^ pRg 

where the quantities with subscript are to be evaluated at the outer edge of the Keplerian 
flow (i.e., at radius r ). 



3.2. Energy Equations for the Keplerian Flow 

In order to get the radial dependence of the flow, we need to solve the energy equation 
to determine the temperature. The gas is heated primarily by the viscous dissipation of 
turbulent energy and Ohmic dissipation (which converts magnetic field energy directly to 
thermal energy), and is cooled by synchrotron and bremsstrahlung emission. Equation ([[]) 
is always valid. By incorporating Ohmic diffusion and viscous dissipation, Equations (^) 
and (0) take the form: 

<H/>V) +pW + pv(V-v) = -L(B-V)B-V (Pth + ^A +2V -(SH , (37) 



^ = V x (v x B) - x (V x B) , (38) 
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respectively, where S {j = (l/2)[v id + v jti - (2/3)5 ij v k)k \ and P th = pR g T/p + P rad is the 
non-magnetic pressure. For a Keplerian flow, S r ^ = —(3/4)0. The radiation pressure, P ra d, 
is given in the Rayleigh- Jeans approximation: 



Prad — —kT 



(39) 



9 v c 

where v m is the frequency below which the radiative emission is highly absorbed, so that 
the optical depth from r to infinity is unity (see Melia 1994). 

Projecting Equation ( |3"TD onto the vector v and Equation (^) onto the vector B, we 
eventually find that 

1 



dt 



P [ 2 v 2 + $ 



-V- 



pv ( -v 2 + $ 



v ■ VP th + 2vV • (Spu) 



v ■ [B x (V x B) 

4:71 



d_ (B 2 
dt \8tt 



1 



737 hH = ^B-[Vx(vxB)] 



2 

CI] 



V • [B x (V x B)] - r)J 2 



An L /J 16n 2 

where J = (c/An)V x B is the current density. Adding the two equations gives 
'B 2 (\ 
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In steady state, the heating term is therefore inferred to be 
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Following the argument by Balbus et al. (1994), we see that the divergence of the viscous 
flux and Ohmic flux is insignificant, and so we can neglect these terms in the following 
discussion. 

In a Keplerian flow, the velocity fluctuations are small compared to the azimuthal 
component of the velocity and it is the correlated fluctuations in the velocity and 
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magnetic field components that produces the anomalous viscosity (Balbus et al. 1994). So 
v = v^e^ + 5v, with 5v much smaller than = \jGMjr. Introducing this into Equation 
([43|) and neglecting the high order terms, we get the steady state heating rate 

1 d fv r Hr < B 2 > 



-v ■ VP t , - 



Hr dr 



An 



_L_d_ 

Hr dr 
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where the correlated fluctuation on the right hand side is given by 

W rih 
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(44) 
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Next, the temperature can be determined by solving the thermal energy conservation 
equation: 

d(pe) 



dt 



-V • (vpe) + T - A - P th V ■ v 



(46) 



where A is the cooling rate, and pe = ankT + 3P ra d is the thermal and radiation energy 
density. In the fully ionized but non-relativistic limit, a = 3, whereas in the relativistic 
electron limit, a = 9/2. Adding Equations ( f£2|) and fl46|) , we have the complete energy 
conservation equation of the system 

'B 2 (\ 



d_ 

Of 



\n H \2 



-r,r - 2upS 2 - A - V- 



pv 



p 2 



-2z/pv • S + — B x (v x B - x B) 

An Ait 
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In steady state, the derivative with respect to time equals zero. So using the mass 
conservation Equation (j3~5|) for a Keplerian flow and Equations (|4j) , fl45|) , (j3l"l) and ([34]) , 
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we have 
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The temperature equation follows from this once we substitute for H, e, P t h etc., whence 
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The 's denote derivatives with respect to r. 



Now, from Equations (|32"D, (|35|) and fl36|), we find that 

V 3 H' 
2T + 2r~ ~ ~H 

0. 



1 v' H' pf 
- + — H h — 

r v r H p 



T 1 fiy/GMv. 



Thus, 



T r 2(3 v (3 p R g r^T v, 
7 3T' pVGMvr 



p 2r 2T 2f3 u f3 p R g r^fT ' 

Therefore eliminating p' in Equation (|49|), we get a differential equation for the 
temperature: 
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where 



9n V c / /i 
2T /^pt^ , / 5 2 




2 r fiv r \2r r 2r 2 ^ i? f 

£ 3 = ^ P»Pp"4> , 327T /^ m \3 

i? 9 T /if r 9n \ c / 

Using the characteristic parameters for Sgr A*, we find that the radiation energy density 
and pressure are always negligible compared to those for the gas within the Keplerian flow. 
So, in Equation fl55j) we will neglect those terms that depend on v m . Thus this equation 
provides the temperature profile throughout the inner region once the outer boundary 
conditions are specified. 



3.3. Calculation of the Spectrum 

The flux density (at earth) produced by the Keplerian portion of the flow is given by 

^ = ^/^ (1-7) ^ dA, (57) 

where D = 8.5 kpc is the distance to the Galactic Center, u is the observed frequency at 
infinity and v' is the frequency measured by a stationary observer in the Schwarzschild 
frame. (For simplicity, we here assume the metric for a non-spinning black hole. A 
more thorough exploration of the parameter values, including the black hole spin, will be 
discussed elsewhere.) The frequency transformations are given by 



i/y/1 - r s /r , (58) 



V = "l-(Vc)cos* ' (59) 
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where v is the frequency measured in the co-moving frame, and 6 is the angle between the 
velocity and the line of sight. Since the radial velocity is always much smaller than v^, 
we ignore this component in the transformation equations. So cos 6 = sin i cos 0, where % is 
the inclination angle of the axis perpendicular to the Keplerian flow, and is the azimuth 
of the emitting element. When the Doppler shift is included, the blue shifted region is 
located primarily near = while the red shifted region is at = ir. The other quantities 
that are necessary for an evaluation of the flux density are the area element 



and the specific intensity 



dA = = cosi r dr d(f) , (60) 

a/1 - r s /r 



I v , = B> v ,{l-e- T ), (61) 



where 




K> = , v , 7-v— - Q B v , (62) 



and the optical depth is 

r=Uds = Kv ^- 1 ~^ Cm \ (63) 

J COSt ^/l _ v 2 /c 2 

where k u is the absorption coefficient. In the case where the optical depth r < 1, Kirchoff's 
law allows us to write 

2 



2H 



I v ,^B' u ,T = e v : V , (64) 

cos^ I 1 — [y^/c) cosu 



where e u = B v k v is the emissivity. The presence of a substantial azimuthal component 
of the magnetic field makes it convenient to calculate the observed flux directly from the 
Extraordinary and Ordinary components of the intensity. The most convenient approach is 
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to select the symmetry axis of the Keplerian flow as the reference direction. The observed 
flux densities in the azimuthal and the reference directions are given by 

F luo = -L J {II, C os 2 0' + 7° sin 2 0') (l - ^) 3/2 dA , (65) 
F 2uo = J [II, sin 2 0' + C cos 2 0') (l - ^) cL4 , (66) 

respectively, where 0' + 7r/2 is the position angle of the magnetic field vector within the 
emitting element that has an azimuth of 0, so that cot0' = cot0 cos?. I*, and 1°, are the 
specific intensities for the Extraordinary and Ordinary waves, respectively. For thermal 
synchrotron radiation, the emissivities are given by (Pacholczyk 1970) 



poo 

„-Bsin0' / N(E)[F(x) + G(x)} dE , (67) 
<,-nm e c z J 
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l e<- Jo 

where N(E) is the electron distribution function at energy E, and 

_ co S fl-y c 

S " l-(v*/c)<xx0' [ ' 

3eB sin 6'E 2 ' 1 ' 

poo 

F(x) = x K 5/3 (z) dz , (71) 

J X 

G(x) = xK 2/3 (x). (72) 

K 5 / 3 and K 2 /s are the corresponding modified Bessel functions. The total flux density 
produced by the Keplerian portion of the flow is the sum of these two. 

Given the temperature profile indicated in Figure 2, a non-negligible fraction of the 
radiation with frequency higher than ~ 3 x 10 11 Hz is self-Comptonized by the hot plasma 
into the X-ray band. We note that for the conditions prevalent in this region (notably 
the temperature of the electrons and the characteristic energy of the seed photons), the 
electron-photon scattering occurs in the Thomson limit. Since in addition the electron 
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thermal energy (~ kT) greatly exceeds the kinetic energy associated with the radial and 
Keplerian motion of the gas, the electrons can be treated as locally isotropic. In turn, 
the upscattered radiation will also be emitted isotropically. Without loss of generality, we 
can therefore calculate the photon number density function n e from the angle-integrated 
intensity (where n t de represents the number density of photons with energies between e and 
e + de and depends primarily on the radius r). Under these conditions, a single electron 
with Lorentz factor 7 emits photons with energies between e s and e s + de s at a rate (e.g., 
Melia & Fatuzzo 1989) 



where \i = cos(#) and fj, s = cos(9 s ) represent the lab-frame (unprimed) propagation angles 
(with respect to the electron's direction of motion) of the incident and scattered photons, 
respectively. Here, we have made use of the transformations e u = e s 7 2 (l — (3fj, s )(l + 0), and 
ei = e s (l — P/j, s )/(1 + (5). The function g(/i',ii f s ) takes into account the angular dependence 
of the cross-section, and is given in terms of the rest-frame (primed) angles by the expression 
#(//, //J = (1 - At' 2 ) (1 - n's) + (1/2) (1 + fi' 2 ) (1 + fj,f). The rest frame and lab frame angles 
are related through the transformation // = (// — /3)/(l — Conservation of energy 
requires that \x = {1 — (e s /e) [1 — P/jl s ]} /P- 

With the upscattered photons emitted isotropically, the total photon production rate 
is then found by integrating over all of the electrons within the emission region, yielding 



where n e is the (radially dependent) number density of electrons and f(E) is the relativistic 
Maxwell-Boltzmann (MB) distribution function normalized to unity. Note that the 
temperature (and therefore the MB distribution function) is also radially dependent. 
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3.4. Results 

Not all the solutions to Equation fl5"5]) are physically acceptable. Equation (|36|) shows 
that for given parameters (3 p and ft v , some boundary values of T , r and tvo will make the 
denominator vanish at a critical radius r crit . For a physically meaningful flow, we therefore 
need r cr u < r,, where r$ is the inner boundary of the Keplerian flow ( = rs in the 
following discussion ). For simplicity, we shall set r cr i t = r^, for which 

This fixes T in terms of r and t> r o- 

A second constraint is provided directly by the current IR and UV observations of 
Sgr A*, which rule out the possible presence of an optically thick disk. Working within 
these constraints, the best fit model for the sub-millimeter bump is that shown in Figures 
1-7. The latest Chandra observation (indicating a luminosity for Sgr A* of no more than 
2.66 x 10 33 ergs s -1 between 2 and 10 keV) restricts the possible range of M near the black 
hole. The reason for this can be understood qualitatively as follows. 

The requirement that the flow remains optically thin in the Keplerian region (i.e., 
within the nominal 5 rs we are considering here) means that naxH < cosi, where 
o~t = 6.65 x 10~ 25 cm 2 is the Thomson scattering cross section. By using Equation fl3"5|) and 
X = 2Hnm p , we therefore have 

M < —47im p rv r cos i / <tt ■ (76) 
But Equations (|75|) and (|36|) show that 

Tr f3 v f3pRg (77) 



/i(GM) 1 /2( r l/2 _ r }/ 2 } 

Combining these two relations, we get 



2II cosi ( r\ 2 ( T \ { Pv \ { P P \ — i 
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Fig. 1. — The spectrum corresponding to the best fit model, whose parameter values are 
indicated within the figure itself. The dashed curve corresponds to the thermal synchrotron 
plus bremsstrahlung component, whereas the dotted curve shows the self-Comptonized 
spectrum. The solid curve is the sum of these. The disk has an inclination angle of 60°. It is 
also necessary to specify the ratio of v r to its free-fall value at r . For this model, this ratio 
is 5.0 x 1CT 4 . The references for the data are given in the text. 
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Fig. 2. — Temperature profile corresponding to the best fit model, shown in Figure 1. 
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Fig. 3. — Magnetic field intensity corresponding to the best fit model, shown in Figure 1. 
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Fig. 4. — Radial velocity profile corresponding to the best fit model, shown in Figure 1. 
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Fig. 5. — Particle number density corresponding to the best fit model, shown in Figure 1. 
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Fig. 6. — Scale height of the gas as a function of r, divided by the Schwarzschild radius, for 
the best fit model, shown in Figure 1. 



-32 - 



Since the factor (r/Vj) 2 / (yjr/ri — 1) has a minimum value of 4 4 /3 3 at r = 16/9 rg, and the 
Keplerian flow must be optically thin everywhere 



M<1.7xl0 21 cosW^-) (V^) gs" 1 . (79) 

Vio 11 k; \o.27 y vo.036; & v ; 

where T m is the maximal temperature. In itself, this value is consistent with the Bondi-Hoyle 
rate inferred from the simulations by Coker & Melia (1997). 

An even stricter upper limit for the accretion rate can be obtained if we assume that the 
sub-millimeter bump in the radio is produced by a hot Keplerian flow through synchrotron 
emission. This plasma will also produce X-radiation via bremsstrahlung emission. The 
bremsstrahlung emissivity is given by 

e b = 6.8 x i -™ n 2 T -^ e -hv/kT ergg cm - 3 g -i Hz -i _ ( g0 ) 

In the Keplerian flow, hu <C kT within the 2 — 10 keV energy range. Thus, with 
exp(—hu/kT) pa 1, we have 

r wkeV d(hv) f ro 

L = / -V^ / dr A7rrHe b = 2.66 x 10 33 ergs s _1 , (81) 

i 2 keV h J r 



and using Equation (|35|), this gives 



r° dr 

M 2 / — < 0.69 x 10 6 (c.g.s.). (82) 



But from Equations ([77]) and fl3"2]) we have (again in c.g.s. units) 
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where xq = ro/Vj. So 
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Fig. 7. — Enlargement of the radio portion of the spectrum for the best fit model shown 
in Figure 1. The dotted curve here corresponds to the first component. The dashed curve 
corresponds to the second component. The values of the parameters are those specified in 
Figure 1. 
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Fig. 8. — Illustration of the effect on the spectrum due to a change in the ratio v r /vff at r . 
In this case, it is 5.0 x 1CT 3 , compared to a value ten times smaller for the best fit model. 
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Fig. 9. — Illustration of the effect on the spectrum due to a change in the value of M, which 
is here five times larger than that of the best fit model shown in Figure 1. 
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We point out here that the range of values giving reasonable fits to the sub-mm data 
(e.g., Fig. 1) falls below this limit, though close to it. Thus, the current high-energy 
observations of Sgr A* do not appear to be in conflict with this model. Indeed, as one 
can see from this figure, the self-Comptonized spectrum fits the Chandra data rather well 
(Baganoff et al. 2000). 

To sample the dependence of our fit on the boundary conditions, we show in Figure 
8 the spectrum resulting from a Keplerian flow forming with a higher radial velocity at 
its outer boundary. Equation ( ff5|) shows the outer boundary temperature inceases by the 
same order. The implied higher viscosity results in a higher temperature over a more 
extended range in radii, producing a much flatter "hump" that does not adequately fit 
the data. This also produces too many X-rays, violating the Chandra limit. However, 
increasing the accretion rate through this Keplerian region will raise the cooling rate so the 
gas temperature decreases. But the increased column density results in a spectral bump 
peaking at higher frequencies than what is observed. In addition, these conditions lead to 
a lower X-ray flux than is observed, though this might still be consistent if the Chandra 
source is not the actual counterpart to Sgr A*. 



4. Concluding Remarks 

Our calculations have shown that under the right conditions, the thermal synchrotron 
emission from a compact, Keplerian flow around Sgr A* can account for the spectral 
(sub-mm) bump. Clearly, the peak frequency of this component is correlated to the gas 
temperature and column density, which characterize several key ingredients, such as the 
accretion rate and the advected specific angular momentum. The latter appears to be 
consistent with the values inferred from the large scale Bondi-Hoyle accretion (Coker & 
Melia 2000), but the required M within the circularized region is significantly lower than 
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that expected further out. This contrast probably indicates the need for a gradual mass 
loss, perhaps in the form of a wind, toward smaller radii. In other words, compared to 
the value of M inferred from the hydrodynamic simulations, the maximum accretion rate 
permitted by the sub-mm data is several orders of magnitude smaller. Yet the specific 
angular momentum with which the gas circularizes is not noticeably different from that 
prescribed by the numerical simulations. 

We note that this model would not pass the Chandra test (of a very low X-ray 
luminosity) if the physical conditions required to produce the sub-mm emission were more 
extreme than those sampled in this paper. Specifically, the sub-mm and X-ray data would 
not be mutually consistent within the context of this model if the accretion rate through 
the Keplerian region was at the Bondi-Hoyle value. 

At this point, we can only speculate about what the overall accretion pattern might 
be like, and we target a future coupling between this analysis and a hydrodynamic 
simulation of the Keplerian flow merged with the larger scale accretion to address the 
issue of self-consistency. We anticipate that when the quasi-spherical infall begins to 
circularize (presumably around 100 — 1000 r$), the turbulent mixing of gas elements with 
high eccentricity readily dissipates the gravitational and kinetic energy densities, raising 
the internal energy of the gas and possibly leading to a rapid evaporation away from the 
rotation plane. This expulsion of a wind may cease when the eccentricity reaches a value 
close to zero, which should occur at a radius ~ 5 — 10 r S - At this point, the flow is very 
nearly Keplerian, and the dynamo process may become active and efficient, along the 
lines we have developed in this paper. The spectrum of Sgr A* longward of 1 — 2 mm is 
presumably generated within the turbulent mixing (i.e., transition) region. 

The first generation of hydrodynamic simulations to address these questions are 
currently underway and we hope to report the results of this analysis in the very near 
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future. 
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